Pauli spin susceptibility of a strongly correlated two-dimensional electron liquid 



A. A. Shashkin*, S. Anissimova, M. R. Sakr^, and S. V. Kravchenko 
Physics Department, Northeastern University, Boston, Massachusetts 02115, U.S.A. 

V. T. Dolgopolov 

Institute of Solid State Physics, Chernogolovka, Moscow District 142432, Russia 

T. M. Klapwijk 

Kavli Institute of Nanoscience, Delft University of Technology, 2628 CJ Delft, The Netherlands 
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Presently, theoretical description of interacting elec- 
tron systems is restricted to two limiting cases: (i) weak 
electron-electron interactions (small ratio of the Coulomb 
and Fermi energies = Ec/Ep <IC 1, high electron 
densities) and (ii) very strong electron-electron inter- 
actions (vg ^ 1, very low electron densities). In the 
first case, conventional Fermi-liquid behavior |lj is es- 
tablished, while in the second case, formation of the 
Wigner crystal is expected (for recent developments, 
see Ref. [2j). Numerous experiments performed in both 
three- (3D) and two-dimensional (2D) electron systems 
at intermediate interaction strengths (1 ^ < 5) have 
not demonstrated any significant change in properties 
compared to the weakly-interacting regime (see, e.g., 
Refs. 0, 0). It was not until recently that qualitative 
deviations from the weakly-interacting Fermi liquid be- 
havior (in particular, the drastic increase of the effec- 
tive electron mass with decreasing electron density) have 
been found in strongly correlated 2D electron systems 
{ts > 10) 1^. However, these findings have been based 
solely on the studies of a kinetic parameter (conductiv- 
ity), which, in general, is not a characteristic of a state 
of matter. 

The 2D electron system in silicon turns out to be a very 
convenient object for studies of the strongly correlated 
regime due to the large interaction strengths (r^ > 10 
can be easily reached) and high homogeneity of the sam- 
ples estimated (from the width of the magnetocapaci- 
tance minima inperpendicular magnetic fields) at about 
4 X 10^ cm~^ In this Letter, we report measure- 
ments of the thermodynamic magnetization and density 
of states in such a low-disordered, strongly correlated 
2D electron system in silicon. We concentrate on the 
metallic regime where conductivity a ^ e'^/h. We have 
found that in this system, the spin susceptibility of band 
electrons (Pauli spin susceptibility) becomes enhanced 
by almost an order of magnitude at low electron den- 
sities, growing critically near a certain critical density 
n-^ ^ 8 X 10^'^ cm~^: behavior that is characteristic in 
the close vicinity of a phase transition. The density 
is coincident within the experimental uncertainty with 



the critical density Uc for the zero-field metal-insulator 
transition (MIT) in our samples. The nature of the low- 
density phase [us < n^) still remains unclear because 
even in the cleanest of currently available samples, it is 
masked by the residual disorder in the electron system. 

Measurements were made in an Oxford dilution refrig- 
erator on low-disordered (lOO)-silicon samples with peak 
electron mobilities of 3 m^/Vs at 0.1 K and oxide thick- 
ness 149 nm. These samples are remarkable by the ab- 
sence of a band tail of localized electrons down to electron 
densities « 1 x 10^^ cm~^, as inferred from the coin- 
cidence of the full spin polarization field obtained from 
parallel-field magnetotransport and from the analysis of 
Shubnikov-de Haas oscillations (the former is influenced 
by possible band tail of localized electrons, while the lat- 
ter is not; for more details, see Refs. (SjiSjiSJ)- This allows 
one to study properties of a clean 2D electron system 
without admixture of local moments |^ 0, ^3 ■ The sec- 
ond advantage of these samples is a very low contact re- 
sistance (in "conventional" silicon samples, high contact 
resistance becomes the main experimental obstacle in the 
low density/low temperature limit). To minimize contact 
resistance, thin gaps in the gate metalization have been 
introduced, which allows for maintaining high electron 
density near the contacts regardless of its value in the 
main part of the sample. 

For measurements of the magnetization, the paral- 
lel magnetic field B was modulated with a small ac 
field i?mod in the range 0.01 - 0.03 T at a frequency 
/ = 0.45 Hz, and the current between the gate and the 
two-dimensional electron system was measured with high 
precision (~ 10""'^^ A) using a current-voltage converter 
and a lock-in amplifier. The imaginary (out-of-phase) 
current component is equal to i = (27r/CBmod/e) dfi/dB, 
where C is the capacitance of the sample and fj, is the 
chemical potential. By applying the Maxwell relation 
dM/dus = —d^/dB, one can obtain the magnetization 
M from the measured i. A similar technique has been 
applied by Prus et al. ^Jl to a 2D electron system in sili- 
con with high level of disorder, in which case the physics 
of local moments has been mainly studied. As discussed 
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FIG. 1; Imaginary current component in the magnetization 
experiment as a function of the electron density in a magnetic 
field of 5 T and T — 0.4 K. Grey area depicts the insulating 
phase. Magnetization versus ris is displayed in the inset. Note 
that the maximum M is coincident within the experimental 
uncertainty with hbtis. 

below, the data analysis and interpretation is not quite 
correct in Ref. in particular, Prus et al. do not dis- 
tinguish between the Pauli spin susceptibility of band 
electrons and the Curie spin susceptibility of local mo- 
ments. 

For measurements of the thermodynamic density of 
states, a similar circuit was used with a distinction that 
the gate voltage was modulated and thus the imaginary 
current component was proportional to the capacitance. 
Thermodynamic density of states dus/dii is related to 
magnetocapacitance via 1/C = I/Cq + l/Ae'^{dns/d^), 
where Co is the geometric capacitance and A is the sam- 
ple area. 

A typical experimental trace of i{ns) in a parallel mag- 
netic field of 5 T is displayed in Fig. ^ The inset 
shows magnetization M(ns) in the metallic phase ob- 
tained by integrating dM/dus — —dfi/dB with the in- 
tegration constant M(cx)) — Bxo, where xo is the Pauli 
spin susceptibility of non-interacting electrons. A nearly 
anti-symmetric jump of i{ns) about zero on the y-axis 
(marked by the black arrow) separates the high- and low- 
density regions in which the signal is positive and neg- 
ative (M(ns) is decreasing and increasing), respectively. 
Such a behavior is expected based on simple considera- 
tions. At low densities, all electrons are spin-polarized in 
a magnetic field, so for the simple case of non-interacting 
2D electrons one expects dji/dB = —^b (at Us — > 0, 
deep in the insulating regime, the capacitance of the sys- 
tem vanishes and, therefore, the measured current ap- 
proaches zero). At higher densities, when the electrons 
start to fill the upper spin subband, M(ns) starts to de- 
crease, and dfj,/dB is determined by the renormalized 
Pauli spin susceptibility x and is expected to decrease 
with Us due to reduction in the strength of electron- 
electron interactions. Finally, in the high-density limit, 
the spin susceptibility approaches its "non-interacting" 
value xoi and d^/dB should approach zero. The onset 
of complete spin polarization — the electron density Up 
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FIG. 2: (a) The experimental dfi/dB as a function of elec- 
tron density in different magnetic fields and T = 0.4 K. The 
curves are vertically shifted for clarity. Grey area depicts the 
insulating phase. Note that the onset of full spin polarization 
in our experiment always takes place in the metallic regime, 
(b) Scaling of the dji/dB curves, normalized by magnetic field 
magnitude, at high electron densities. The dashed line rep- 
resents the "master curve". Spin susceptibility obtained by 
integrating the master curve (dashed line) and raw data at 
B — 1.5, 3, and 6 T is displayed in the inset. 



at which the electrons start to fill the upper spin sub- 
band — is given by the condition dji/dB = {M{ns) 
reaches a maximum), as indicated by the black arrow in 
the figure. It is important that over the range of magnetic 
fields used in the experiment (1.5-7 tesla), the maximum 
M coincides within the experimental uncertainty with 
fiBns thus confirming that all the electrons are indeed 
spin-polarized below Hp. Note however that the absolute 
value of dfi/dB at < Uc is reduced in the experiment. 
We attribute this to smearing of the minimum in i{ng) 
caused by possible influence of the residual disorder in 
the electron system, which is crucial in and just above 
the insulating phase, in contrast to the clean metallic 
regime we focus on here. Another reason for the reduc- 
tion in dfi/dB is the electron-electron interactions (due 
to, e.g., the enhanced effective mass). 

In Fig.[2Ia), we show a set of curves for the experimen- 
tal dfi/dB versus electron density in different magnetic 
fields. Experimental results in the range of magnetic 
fields studied do not depend, within the experimental 
noise, on temperature below 0.6 K (down to 0.15 K which 
was the lowest temperature achieved in this experiment). 
The onset of full spin polarization shifts to higher electron 
densities with increasing magnetic field. Grey area de- 
picts the insulating phase, which expands somewhat with 
B (for more on this, see Ref. Note that the range of 

magnetic fields used in our experiment is restricted from 
below by the condition that dfi/dB crosses zero in the 
metallic regime. In Fig.|21[b), we show how these curves, 
normalized by magnetic field, collapse in the partially- 
polarized regime onto a single "master curve" . The exis- 
tence of such scaling verifies proportionality of the mag- 
netization to B, confirming that we deal with Pauli spin 
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FIG. 3: (a) Magnetocapacitance versus electron density for 
different magnetic fields, (b) Deviation of the C{na) depen- 
dences for different magnetic fields from the B — reference 
curve. The traces are vertically shifted for clarity. The onset 
of full spin polarization is indicated by arrows. 



susceptibility of band electrons, and establishes a com- 
mon zero level for the experimental traces. Integration 
of the master curve over Us yields the spin susceptibil- 
ity X = M/B, as shown in the inset to Fig. |2Ib). Also 
shown is the spin susceptibility obtained by integration 
of raw curves at B — 1.5, 3, and 6 tesla, which, within 
the experimental error, yield the same dependence. 

This method of measuring the spin susceptibility, being 
the most direct, suffers, however, from possible influence 
of the unknown diamagnetic contribution to the mea- 
sured dji/dB, which arises from the finite width of the 2D 
electron layer 13]. To verify that this influence is negligi- 
ble in our samples, we employ another two independent 
methods to determine x- The second method is based on 
marking the electron density rip at which dfi/dB = and 
which corresponds to the onset of complete spin polar- 
ization, as mentioned above. The so-determined polar- 
ization density np{B) can be easily converted into xi^^s) 
via X = ^^B'r^p/B. Note that in contrast to the value 
of dfi/dB, the polarization density Up is practically not 
affected by possible influence of the diamagnetic shift. 

The third method for measuring Up and x, insensitive 
to the diamagnetic shift, relies on analyzing the magne- 
tocapacitance, C. Experimental traces C{ns) are shown 
in Fig. Ola) for different magnetic fields. As the mag- 
netic field is increased, a step-like feature emerges on 
the C(rts) curves and shifts to higher electron densities. 
This feature corresponds to the thermodynamic density 
of states abruptly changing when the electrons' spins be- 
come completely polarized. To see the step-like feature 
more clearly, in Fig.|2Ib) we subtract the C{ns) curves for 
different magnetic fields from the reference B = curve. 
The fact that the jumps in C (as well as in dfi/dB) are 
washed out much stronger than it can be expected from 
possible inhomogeneities in the electron density distribu- 
tion (about 4 X 10^ cm~^ 0) points to the importance of 
electron-electron interactions. Since the effects of inter- 
actions are different in the fully- and partially-polarized 
regimes, it is natural to mark the onset of full spin po- 



FIG. 4: Dependence of the Pauli spin susceptibility on elec- 
tron density obtained by all three methods described in text: 
integral of the master curve (dashed line), dfi/dB — (cir- 
cles), and magnetocapacitance (squares). The dotted line is 
a guide to the eye. Also shown by a solid line is the transport 
data of Ref. Inset: polarization field as a function of the 
electron density determined from the magnetization (circles) 
and magnetocapacitance (squares) data. The symbol size for 
the magnetization data reflects the experimental uncertainty, 
and the error bars for the magnetocapacitance data extend to 
the middle of the jump in C. The data for Be are described 
by a linear fit which extrapolates to a density close to the 
critical density ric for the B = MIT. 



larization at the beginning of the interaction-broadened 
jump, as indicated by arrows in the figure. In case the 
residual disorder does contribute to the jump broadening, 
we extend error bars to the middle of the jump, which 
yields an upper boundary for the onset of full spin polar- 
ization. 

In Fig. ^ we show the summary of the results for the 
Pauli spin susceptibility as a function of n^., obtained 
using all three methods described above. The excel- 
lent agreement between the results obtained by all of the 
methods establishes that a possible infiuence of the dia- 
magnetic shift is negligible and, therefore, the valid- 
ity of the data including those at the lowest electron den- 
sities is justified. There is also good agreement between 
these results and the data obtained by the transport ex- 
periments of Ref. . This adds credibility to the trans- 
port data and confirms that full spin polarization occurs 
at the field B^, however, we note again that evidence for 
the phase transition can only be obtained from thermo- 
dynamic measurements. The magnetization data extend 
to lower densities than the transport data, and larger 
values of x are reached, exceeding the "non-interacting" 
value xo by almost an order of magnitude. The Pauli 
spin susceptibility behaves critically close to the critical 
density ric for the B = metal- insulator transition (15| : 
X oz ris/ (ris — n^). This is in favor of the occurrence of a 
spontaneous spin polarization (either Wigner crystal 
or ferromagnetic liquid) at low rij, although in currently 
available samples, the formation of the band tail of lo- 
calized electrons at rig < ric conceals the origin of the 
low-density phase. In other words, so far, one can only 
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reach an incipient transition to a new phase. 

The dependence Bdng), determined from the magne- 
tization and magnetocapacitance data, is represented in 
the inset to Fig. 0] The two data sets coincide and are 
described well by a common linear fit which extrapolates 
to a density close to Uc- We emphasize that in the 
low-field limit {B < 1.5 tesla), the jump in dji/dB shifts 
to the insulating regime, which does not allow us to ap- 
proach closer vicinity of n^: based on the data obtained 
in the regime of strong localization, one would not be 
able to make conclusions concerning properties of a clean 
metallic electron system which we are interested in here. 
Clearly, the fact that the linear Bc{ns) dependence per- 
sists down to the lowest electron densities achieved in the 
experiment confirms that we always deal with the clean 
metallic regime. 

Finally, we would like to clarify the principal differ- 
ence between our results and those of Ref. [1^. In the 
sample used by Prus et al, the critical density for the 
B = MIT was considerably higher than in our sam- 
ples caused by high level of disorder, and the band tail 
of localized electrons was present at all electron densi- 
ties As a result, the crucial region of low electron 
densities, in which the critical behavior of the Pauli spin 
susceptibility occurs, falls within the insulating regime 
where the physics of local moments dominates 0, Ifllllflll . 
Indeed, Prus et al. have found sub-linear M{B) depen- 



dence characteristic of local moments, and the extracted 
spin susceptibility in their sample has a Curie tempera- 
ture dependence |^ . This is the case even at high electron 
densities, where metallic behavior might be expected in- 
stead. Such effects are absent in our samples: the spin 
susceptibility (in the partially-polarized system) is inde- 
pendent of the magnetic field and temperature, confirm- 
ing that we deal with Pauli spin susceptibility of band 
electrons. 

In summary, the Pauli spin susceptibility has been de- 
termined by measurements of the thermodynamic mag- 
netization and density of states in a low-disordered, 
strongly correlated 2D electron system in silicon. It is 
found to behave critically near the zero-field MIT, which 
is characteristic of the existence of a phase transition. 
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